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1 Introduction 



The study of the quantum properties of grand unified theories (GUTs) in the presence of 
a strong curvature is quite an interesting issue, owing to the different apphcations that it 
can have in a number of situations. The resuhs of this study are important for a detailed 
knowledge of the early universe, in particular for an accurate discussion of the known models 
of inflationary universe (see [||] for a review) and as a guide in the construction of new models 
of this kind. Moreover, such a study is relevant for a better understanding of quantum effects 
in the vicinity of a black hole — in particular, for instance, for the calculation of quantum 
corrections to the black- hole entropy (a recent discussion can be found in 0). Furthermore, 
such considerations are of fundamental importance for the estimation of the back reaction 
effect a quantum field has on the geometry of spacetime (for an earlier discussion, in the free 
matter case, see |0, ^). 

The renormalization group (RG, see ||^ for an introduction) turns out to be very useful 
in the discussion of GUTs in curved space (see, for example, 0). The first investigations 
on these topics, which included the construction of the RG for GUTs in curved spacetime 
(see [0, H, and for a complete list of references), have been followed by a lot of activity, 
where the subject is considered under quite different points of view. Among the different 
interesting phenomena which are specific of this theory we can count curvature-induced 
asymptotic freedom 0, asymptotic conformal invariance both in asymptotically- free ^ ^ 
and in finite GUTs applications of phase transitions of Coleman- Weinberg type in 
inflationary universes curvature-induced phase transitions [|IT|, and so on. 

In the present paper we study the renormalization group properties of GUTs in curved 
spacetime. We start from asymptotically free GUTs in curved spacetime and write the whole 
system of RG equations for all the gravitational couplings. Their behaviour is not asymptot- 
ically free, of course. Concentrating mainly on the gravitational coupling constant G, we give 
its general running form in asymptotically free GUTs and provide some explicit examples 
for the gauge groups SU{2), Eq (Sect. 2). We also discuss the running couplings for finite 
GUTs in curved spacetime. The running gravitational coupling is calculated explicitly and it 
is shown that quantum corrections to G have an exponential form, unlike in asymptotically 
free GUTs, where they behave power- wise (Sect. 3). 

Sect. 4 is devoted to the use of the running gravitational constant for calculations of 
radiative corrections to the Newton potential. In Sect. 5, in order to study how QG effects 
may change the qualitative picture obtained in the previous section, we discuss the effective 
theory of conformal gravity by Antoniadis and Mottola. This theory aims at the description 
of IR quantum gravity. The running gravitational constant in such a theory looks qualita- 
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tively similar to that in asymptotically free GUTs, but different from the one in the Einstein 
theory, where quantum corrections to the Newtonian potential have also been calculated 
recently P3| . 

Sect. 6 is devoted to the construction of a totally asymptotically free theory of matter 
with i?^-gravity. Considering an 0{N) gauge theory with one multiplet of scalars and two 
multiplets of spinors as matter, and making use of the very interesting technique of intro- 
ducing potentials in RG-coupling space (similar to a c-function), we explicitly construct the 
regimes where the total matter-QG theory is asymptotically free (such study is carried out 
numerically). Then, in Sect. 7 we compare the behaviour of the running gravitational cou- 
pling in this matter-QG system with the running of G in the same theory with QG being 
classical. Sect. 8 is devoted to the study of quantum corrected gravitational equations in the 
strong curvature regime. The non-singular De Sitter solution of these equations with GUT 
corrections is discussed, as well as a solution of wormhole type. In the concluding section 
we summarize our results and outline some possible extension of the same. 

2 Asymptotically-free GUTs in curved spacetime and 
the gravitational coupling constant 

Our considerations start from a specific GUT in curved spacetime, given by the following 
Lagrangian (a multiplicatively renormalizable one, H) 



1(^ 



'.ext = aiR^ + a2Cl^^. + a3G + a4nR + A- -^R. (2.1) 



With an appropriate gauge group, the theory defined by the Lagrangian (|2.1| ) contains gauge 
fields A^, scalars ip and spinors ip, in some representation of the given gauge group. As usual, 
the Lagrangian of the external fields must be added to in order to obtain a theory which 
is multiplicatively renormalizable in curved spacetime [^]. In ( p.l|) the cosmological constant 
has been chosen in a specific form which will be convenient in order to facilitate the discussion 
of the cosmological applications below. 

The detailed consideration of the renormalization structure and the RG equations for an 



asymptotically free GUT of the form (|2yi| ), based on the gauge groups SU(5), SU(2), 0(N), 
Eg (and some others) can be found in the book where relevant references are listed 
too. Now, since the theory is multiplicatively renormalizable, the effective Lagrangian that 
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corresponds to the classical theory (|2.1|) fulfills the standard RG. We shall assume that the 
background fields vary slowly with respect to the effective mass of the theory and, therefore, 
the derivative expansion technique can be used in order to obtain the effective action of the 
theory. Then, we can restrict our considerations to the terms with only two derivatives with 
respect to the scalar fields and four derivatives with respect to the purely gravitational terms, 
and it turns out that the structure of the effective Lagrangian just mimics the structure of 
the classical Lagrangian 

The RG equations satisfied by the effective Lagrangian are: 

( d d d \ 

T^^eff = ( + Pi-Q^ - '^i^i'o^.j ^effilJ', K, </>i) = 0, (2.2) 

where /i is the mass parameter, Aj = ((7^, /i^, /, m^, ^, A, G, ai, 02, 03, 04) is the set of all 
coupling constants, the jSi are the corresponding beta functions, and (pi = [A^, (p, ip) are 
the fields. Note that the dependence on the external gravitational field g^j^^ is not explicitly 
shown in L^q^. 

The solution of Eq. ( p.2|) by the method of the characteristics gives (for a similar discus- 
sion in the case of the RG improved Lagrangian in curved space, see [0): 

Le//(/i, Ai, (pi) = Leff{ne\ Xi{t), (piit)), (2.3) 

where 

^ = A(A.(t)), A.(0) = A„ 

^ = -7^m.{t), 0.(0) = 0.. (2.4) 

The physical meaning of (( |2.3| )) and ((|2.4|)) is that the effective Lagrangian L^fj (called 
sometimes the Wilsonian effective action |^) is found provided its functional form at some 
value of t is known (usually the classical Lagrangian serves as boundary condition at t = 0). 
We will come back to the discussion of L^ff later on and, for the moment, we shall concentrate 
on the scaling dependence of the coupling constants, Eq. ( p.4|) . Note that only using the RG- 
improvement procedure can one also get the non-local effective action, which was discussed 
in Refs. ^ by direct one-loop calculation. 

We consider a typical asymptotically free GUT in curved spacetime. For studying the 
scaling dependence of the coupling constants, the RG parameter will be chosen to be t = 
ln(/i//io), as usually, where /i and /iq are two different mass scales (see for a rigorous 
discussion of the RG in curved spacetime). The running coupling constants corresponding 
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to asymptotically free interaction couplings of the theory have the form 



9\t) = / 1 + 



\ (4vr)2 J ' 
h\t) = K,g\t), fit) = K,g\t) 



9\0) = g^ 



(2.5) 



where Ki and K2 are numerical constants defined by the specific features of the theory under 
consideration (see [|l^-||l5l for explicit examples of such GUTs in flat space, and ^ for a 
review). As one can see, asymptotic freedom (g'^it) — 0, t — 00) is realized [l^, for the 
gauge coupling and for the running Yukawa and scalar couplings. 

The study of this kind of asymptotically free GUT in curved spacetime was been started 
in Refs. [0, §] (for a review see 0), where the formulation in curved space was also developed. 
Using the results in those works it is easy to show that in such theories (for simplicity, with 
only one massive scalar multiplet), one obtains 



1 

6 



B'gH 



m 



1 + 



B'gH 



[2.6) 



being ,^(0) = ^, m^(0) = m^, and where for the different models the constant b can be either 
positive [0 H, P], negative 0, 0], or zero 0. Notice that the constant one-loop running 
coupling ^(t) = ^ (i.e. 6 = 0) usually corresponds to supersymmetric asymptotically free 
GUTs, and also that classical scaling dimensions are not included in the RG equations for 
couplings with mass dimension, as usually happens in the RG improvement procedure. 
Turning now to the gravitational coupling constants @] (see also [0, ||) 



dai{t) 

dt 
da2{t) 

dt 
da^lt) 

dt 
dA{t) 

dt 



(47r)2 
1 



m- 



6 



120(47r)2 
1 



12A^, 



d 



360(47r)^ 
m\t)Ns 
2(47r)2 ' 

m'^{t)Ns 



- {Ns + llNf + Q2Na) 



dt IQnGit) 



(4- 



(2.7) 



where A'^,, Nf and Na are, respectively, the number of real scalars, Dirac spinors and vectors 
and we work in euclidean region. In what follows we shall consider dynamical spacetimes 
with constant curvature (as the De Sitter space), so that the term in Di? will not appear. 
As one can see, the behavior of the gravitational couplings a2{t) and 03 (t) is given by 
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02,3 (^) — ^2,3 + 02,3 1, and it is exactly the same as in the free matter theory with the 
same field content. It is not infiuenced by the interaction effects. 

The most interesting quantity for us will be the gravitational running coupling constant 
G. As one can easily find from ( (|2.7|) ), it has the form 



G{t) 



Go 1 



Go 1 + 



167rAr,Gom^(e- 1/6) 
5V(2& + 1) 

IQ-KN.Gom^^i - 1/6) 
SV(26 + 1) 




2b+l 



- 1 



2b+l 



Gn 



1 + — Gom2(^-l/6)t 

TT 



where t = ln(yu/yUo). As one can see, the matter quantum corrections to the gravitational 
coupling constant in GUTs are larger than in the free matter theory if 6 > 0. 

Let us now give some examples. First of all, we consider the asymptotically free SU(2) 
model of Ref. [|T^ with a scalar triplet and two spinor triplets. Then, one can show ||^ that 
Ng = 3, = 10/3, 6 ~ 2. In this model (with the standard choice ~ 0.41) we get for the 
first non-trivial correction to the classical G: 



G{t)-Go 1 + 0.9549 Gom2(^- 1/6) t 



(2.9) 



As a second interesting example we will consider the asymptotically free Eg GUT in 
curved spacetime [IS]. This theory contains a 78-plet of real scalars and two 27-plets of 
charged scalar, M and A^. Assuming that only one mass in the real scalar multiplet is different 
from zero, and choosing the initial values for the charged scalars to be C,m = ^at = 1/6 we 



obtain again a G{t) of the same form as 
52 = 32, g^ = 0.41, we get 



8^. With the following parameters 



1 

6 



0.97(e^ - 1/6) I 
■ B^gH' 



1 + 



1/2 



B'gH 
(47r)2 



1/2 



0.97ml 1 + 



{Any 



(2.10) 



Then, we obtain 



G{t) 



Go a 



16nNsGoml0.97\^^ - 1/6) 
5V(2fe+ 1) 



B'gH 



26+1 



Go 
Go 



l + ^0.972GomJ(e^-l/6)t 
1 + 0.8985 Gom2(e<^- 1/6) t 



(2.11) 
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We thus see that we are able to obtain the general form of the scaling gravitational cou- 
pling constant in an asymptotically free GUT in curved spacetime. As we observe from 
these explicit examples, there exist asymptotically free GUTs in which the matter quantum 
corrections to the gravitational coupling constants at high energies are bigger than these 
corrections in the free matter theory. The application of the above results will be discussed 
later. The gravitational constants for other asymptotically free GUTs can be obtained in a 
similar way. 



3 Finite GUTs in curved spacetime and the gravita- 
tional coupling constant 

A very interesting class of GUTs in curved spacetime is given by finite GUTs (see, for 
instance, [ll9|-|^). Finite GUTs in curved spacetime were first considered in Ref. [|lOl (for 
a review and a complete list of references, see again 0). Of course, although finite for the 
interaction couplings in flat space, these theories are not finite in curved spacetime, due to 
vacuum, mass and Rip'^-tjpe divergences. Recently, quite a spectacular development has 
emerged which concerns N = 2 supersymmetric models (in particular, finite) with matter 
multiplets, through the definition of the exact spectrum [^. For such theories we have (in 
the following we consider one-loop finite supersymmetric or non-supersymmetric theories) 



(3.1) 



and 



rn^it) = exp(C(7^t), 



(3.2) 



where the constants ki and k,2 depend on the specific features of the theory, and where C can 
be positive, negative or zero [l^. Solving the RG equations for the gravitational coupling 
constant in such theory, we find 

-1 



G{t) 



Gn 



Gn 



1 - 



1 + 



N,Gom\^ - 1/6) 
NsGom\^-l/6) 



_2CgH 



1 



(3.3) 



Hence, as one can easily see, GUTs with a positive C are the ones which give the biggest 
contribution to the gravitational coupling constant, among all finite GUTs. 

As an example, one can consider the SU(2) finite model of the first of Refs. [^, with a 
SU(4) global invariance and the scalar taken in the adjoint representation of SU(2). (Notice 
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that there is = 2 supersymmetry here, in one of the regimes of finiteness.) In this case, 
we get that [jlO| A"^ = 18 and C = 2A/{A'nf. As a resuh, 

Gom\i - 1/6) 



G{t) = Go 



1 + 67r 



(e3.^*M^ - l) 



(3.4) 



As one can see, we have power corrections — of the form ~ (^^^ — to the gravitational 
couphng constant. In a similar way one can also obtain the running gravitational coupling 
corresponding to other finite GUTs. 



4 Quantum GUTs corrections to the Newtonian po- 
tential 

As an application of the results of the previous discussion, we will consider in this section 
the quantum corrections to the gravitational potential. From the discussion in Sects. 2 and 
3 we know that the typical behavior of the gravitational coupling constant, when taking into 
account the quantum corrections is 

Git) ~ Go {l + Go [(1 + St)2^+i - l] } , (4.1) 

in asymptotically free GUTs, and 



G{t) Gc 



1 + Go ( e^* - 1 



(4.2) 



in finite GUTs. Of course, such corrections are too small to be measured explicitly, although 
they are significantly larger than in free matter theories (due to running). In the above 
scaling relations, t = In^fi/^o). 

It was suggested in Refs. |2^, ^ that in the running gravitational and cosmological con- 
stants /i//io ought to be replaced with ro/r, i.e. that one should change the mass scalar ratio 
by the inverse ratio between distances. The support for this argument comes from quantum 
electrodynamics, where the well-known electrostatic potential with quantum corrections can 
be alternatively obtained from the classical potential by interchanging the classical electric 
charge with the running one, with t = ln(ro/r). Similarly, one can estimate now the GUTs 
quantum corrections to the gravitational potential. Starting from the classical Newtonian 
potential 

Gmim2 , . 

V{r) = , (4.3) 

r 

the classical gravitational constant in ([4.3|) is to be changed with the running gravitational 



coupling constant (|0 



8 



Gomim2 



As a result, we obtain, 
V{r) '. 

in asymptotically free GUTs, and 

Gomim2 
V{r) ~ 



1 + ciGom^ - ^) In 



1 



Cg^ 



QJ \r J 



(4.4) 



(4.5) 



in finite GUTs. Here Go is the initial value of G (the value at distance vq), and Ci and C2 
are some constants. 



In Ref. |]2^ the gravitational potential has been calculated in the frame of quantum Ein- 
steinian gravity considered as an effective theory (this is unavoidable, owing to its well-known 
non-renormalizability ||26|). It was found that in such theory the leading quantum corrections 
are proportional to Gq/t^. As we see, in an asymptotically free GUT the leading correction 
is of logarithmic form ~ ln(ro/r) while in a finite GUT it is power-like ~ (ro/r)*^^^. Hence, 
from the examples of Sects. 2 and 3 one understands that the GUT quantum corrections to 
the Newtonian potential can be more important than the corresponding corrections in the 
effective theory of Einstein's quantum gravity. We now turn to the study of the gravitational 
constant in some model of quantum gravity. 



5 The gravitational constant in the effective theory for 
the conformal factor 

In order to see how quantum gravity (QG) effects change the results of the above discussion, 
we have to consider some multiplicatively renormalizable QG (with matter). For example, 
in multiplicatively renormalizable pure higher derivative gravity (for a review, see 0) the 
running gravitational coupling has been used to construct the QG corrected Newtonian 
potential in Ref. [Q. It has been suggested there that the QG corrections in the Newtonian 
potential (which are also of logarithmic form) might help to solve the dark matter problem 
and some other cosmological problems p7[] . 

In this section we are going to discuss the running gravitational coupling constant in the 
effective theory for the conformal factor ||2^ , which presumably describes the infrared phase 



of QG. The construction of the effective theory for the conformal factor proceeds as follows 
p8| . One starts from the conformally invariant matter theory (the free theory, for simplicity) 



in curved spacetime. The standard expression for the conformal anomaly is known (see p9| , 
and for reviews 0, 13]), and it can be integrated on a conformally-flat background (in the 
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conformal parametrization g^y = e^'^Qniy), in order to construct the anomaly- induced action. 
Adding to this action the classical Einsteinian action in the conformal parametrization, we 
get the effective theory for the conformal factor. In the notations of Ref. and for the 
flat background g^i, = r]^^, the resulting Lagrangian is 

' - - C \2a{d^a)^na + a\d^a)^] + 7e'"'^(9^a)2 - Ag^-, (5.1) 



where 7 = 3/(87rG'), A = A/(87rG), ( = b + 2b' + 3b" and gV(47r)2 = (- 2b', being b,b' 
and b" the well-known coefficients of the conformal anomaly (for a recent discussion, see 
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= b (Cl^^^ + ^Di?) + b'G + b"uR. (5.2) 

has been interpreted as the 4-dimensional central charge and is the anomalous scaling 
dimension for a. 

Near the infrared (IR) stable fixed point C = 0, the theory ( p.l[ ) has been argued to 
describe the IR phase of QG |^8|]. Near the IR fixed point C = 0? the inverse of the running 
gravitational coupling constant in the IR limit {t —>■ —00) has been calculated in Ref. pT| : 



7(t) ~ (-t)"^/^exp 



t 2 - 2a + 



2 



(5.3) 



where Ql = Q^(C = 0). As a result, we obtain the Newtonian potential with account to the 
quantum corrections in such a model, as [pTl 



Gomim2 ( ^ o? 

logf ""r^ 



V{t) = -^^^ 1 _ in^ . (5.4) 



As we can see here, essential leading-log corrections to the Newtonian potential appear in 
this model, what is different from what happens in the case of Einstein gravity. 



6 Asymptotic freedom in GUTs interacting with higher 
derivative quantum gravity 

After the above discussion on the running gravitational coupling constant, mainly for GUTs 
in curved spacetime, our purpose will be to look at its quantum gravitational corrections in 
frameworks of asymptotically free theories. To be specific, we start from a model with the 
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following Lagrangian (here we use notations slightly different from those in Sect. 2) 
C 



^ ^ C^iivafi ^ jr)2 ^ (~ia /^afiu 



1 



R 



A 



(6.1) 



2'^' levrG SttG 



We will not discuss the problem of unitarity in this theory, which is still open and may 
be solved, perhaps, only non-perturbatively. We rather consider (|6.1| ) as an effective theory 
for some unknown consistent QG. This model includes in its gravitational sector higher 
derivative gravity (with somehow different notations for the gravitational couplings) and in 
the matter sector an 0{N) gauge theory with scalars in its fundamental representation plus 
rii spinor mutliplets in the adjoint representation and n2 spinor mutliplets in a fundamental 
representation. It is known ||I5|, ^ that for some values of N,ni,n2 such a fiat GUT 
provides asymptotic freedom for all its coupling constants. Moreover, even taking into 



account quantum gravity, one can show that asymptotic freedom may survive (see |Q) but 
under tighter restrictions on the contents of the theory. When studying the asymptotically 
free regime we look at the region h"^ ^ g"^ [^] (so that Yukawa coupling corrections may be 
dropped). One can show that such a description is consistent under radiative corrections. It 
is interesting to remark that in the low-energy limit the theory (|6.1|) leads to the standard 



Einstein action [10| — similarly to what happens in the pregeometry program [46|. 

The one-loop beta functions for the theory (|6.1| ) have been calculated in |^3|, |] . We will 
be interested in the study of that system of equations in the t —>■ oo (or high-energy) limit. 



where it is convenient to make the variable changes |33 



/ 

dr 



f 



-.2 ' 



Kg'^{t)dt, 



K = 



{Any 



Further, we will take advantage of the property that, in the limit considered. 



m 



(6.2) 



62 

2 ■ 



9'{t) a 

First, we look at the case of massless couplings only. In these conditions, the total system 
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of RG equations turns into 



where 



dr 

duj 



dr 
d^ 

dr 



dl 
dr 



K 

a" 
62 



2^2 



-1,'g 



4a 



10 



+ (5 + 



12 2 



A^ + 2- 3(A^-1^ 



-/ 



-9 



2 

^ 3.2 o 10 9 ^2 5 ^ 



9 



2cu 



1 

3uj 



+ f 15 + 



4cj2 



4 

27^2^ 



6. 1 



1 



(47r)2' 

^(798 + 6iV2-5iV) + ^ 
^(22iV-45)-^(ni(iV-2)+n2 

O 6 



(6.3) 



(6.4) 



It is evident that asymptotic freedom (AF) in the original couphngs is determined by the 
existence of stable fixed points for this new system. Note that the RG equations ( |6.3D are 
also very useful for explicit discussions of different forms of effective potential in quantum 
matter- i?2-gravity theories |ll], ^ 



6.1 Without gravity 

To make our discussion easier, we start by considering the same situation when gravity is 
switched off. In these conditions, 



dr 
dl 

dr 



1,2 2 
iV + 8.2 



9 



(6.5) 



/2 + (^2_3(^_l))J+ 



4 



This system will be examined using the methods developed in ||3^, where some potentials in 
RG coupling-space were introduced, so that their stability properties yield the fixed points 
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of the original system. One introduces a renormalization group potential u{g , /) such that 



2 



du dg 

dg"^ dr 

du df 

df ^ d^' 



(6.6) 



With the sign convention here used, the existence of stable fixed points amounts to the 
presence of some sort of maximum for u (if the signs in ( |6.6|) were reversed, we should say 
minimum instead of maximum). Up to an arbitrary constant, this potential reads 

First, we obtain the values of the fixed points, which correspond to the critical points of u 
in ((yf^, /)-space. Trivially, g^ = 0, but there are real solutions for / only when 

A = {b'^ - 3{N - 1))2 - 3{N + 8)(iV - 1) > 0. 

If A > 0, the critical values of / are 



2/3(A^ + 8) 

Therefore, the Hessian matrix at these points is 



/. . -'"'-^'f-.^r^^ . (6.8) 



-62 

±Va 



(6.9) 



Obviously, its eigenvalues are 6^ and ±vA. We need 6^ > in order to ensure AF in g'^. 
Thus, u can have an extreme only when picking the minus sign, which corresponds to / = /2, 
and then that extreme is a local maximum. In such a set-up AF for / takes place. By the 
same argument we conclude that / = /i corresponds to a saddle point. The hypothetical 
case A = is rather exceptional; there is only one solution for / and the outcome is one 
vanishing eigenvalue (i.e. the Hessian no longer has maximal rank). Actually, while there 
is a maximum along the g^ axis, the potential shows just an inflexion point along the / 
direction, typical of the cubic dependence of u in this variable. 

When considering only the /-dependent part of u — say uj — we realize that it has no 
global maximum or minimum, because cubic curves are unbounded from above or below. 
However, if A > 0, it has a local minimum (at / = /i) and a local maximum (at / = /2). 
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If A = 0, it has an inflexion point at the only zero of (3f and no local extremes. On the 
contrary, when A < 0, is a monotonic curve without extremes or inflexion points. In the 
end, it is the presence of the local maximum / = /2 that accounts for AF. 

As we have just seen, the AF scenario is determined by 6^ > and A > 0. Whether 
these constraints can be met or not depends, of course, on the particular values of and 
of b which, in turn, depends on A^, rii and n2. Assuming N > 2, and bearing in mind that 
77-1, n2 are integers larger than or equal to zero, from ( |6.4| ) one readily finds the condition 

ni < 2. (6.10) 

Then, after numerical examination of A for 1 < A^ < 10, < ni < 2 and < ^2 < 20, we 
find the following combinations yielding AF solution: 



A^ 


ni 


n2 


7 





13 




1 


8 




2 


3 


8 





15,16 




1 


9,10 




2 


3,4 


9 





17,18,19 




1 


10,11,12 




2 


3, 4,5 


10 





19,20 




1 


11,12,13 




2 


3, 4,5 



This constitutes an example of (discrete) numerical boundary separating the region where 
the theory shows a given behaviour for r ^ cxd. 

Let's now see, from another viewpoint, the meaning of our constraints > and A > 
when looking at the couplings themselves. In the present case, we can integrate the RG 
equations ( |6.5| ). For g"^ one has 

g\r) = g\0)e-''^ (6.11) 

(we are taking initial conditions so that r = when t = 0). It is plain — as we have said — 
that 6^ > leads to AF in g'^, while 6^ < entails boundless increase of (7 as r grows, making 
perturbation theory no longer valid. With regard to /, there are three cases depending on 
A: 
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1. A < (no fixed points) 
fir 



2{N + 



v/Za tan ( ^^^^] - {I? - 3(iV - 1)) 



(6.12) 



/ is periodic in r, having regularly spaced singularities — analogous to Landau poles. 
2. A = (one double fixed point) 

(62-3(Ar-l)) 1 



There is stability in the sense that /(r 
have a pole of this coupling at r = 0. 

3. A > (two single fixed points) 



2/3(iV + 8) 2/3(iV + 8)T 
(62-3(iV-l)) 



(6.13) 



2/3(A^ + 8) 



as r ^ oo, but we still 



f{r) 



{b^ - 3{N - 1) - VA)e-^^ - (6^ - 3{N - 1) + VA) 



2/3{N + 8){l-e-^^) 
It is clear that now a finite limit exists when r goes to infinity. In fact, 

(62-3(iV-l) + v^) 



(6.14) 



fir) 

T 



2/3(iV + 



/2 



(6.15) 



oo, 



as should be expected, because we already knew that /2 is the stable fixed point. Now, 
there can be no doubt that the previous criteria are right and this is indeed the AF 
region. 



6.2 With gravity 

Next, we will have to deal with Eqs. (|6.3|) , which offer more serious difficulties than the 
system (|6.5|) . Since not all the crossed derivatives coincide, we cannot just integrate and 
find a potential which is a function of all the variables. Instead, we may handle individual 



dUg2 




du^ 


dr ' 
du 


§^ 




dUf 


dr^ 
dl 


df 


dr 



(6.16) 



15 



Some possible solutions are 
62 



Ug2 



Ui 




9 ^4^A^3_^^^ 

16to 6iJ 6uj 

6 1 

Zoo uj 2uj 



(6.17) 



Next, we obtain their critical points (which are zeros of the beta functions ( |6.3| )) by numerical 
methods. Note e.g. the difference between the /-coefficients in without QG and in the 
above expression. Is is not difficult to realize that the QG contributions may get to change 
the balance which makes A negative or positive, bringing about modifications in the limits 
of the AF region, as we shall see below. 

Once we have the numerical values of all these fixed points, we classify them according 
to the criterion of the eigenvalues of the matrix 



/ 





d'^Ug2 


d'^Ug2 


d'^Ug2 




dg'^duj 


dg^d^ 


dg^df 






duodS, 


dujdf 


d^dg^ 


d^duj 
d\f 


de 


dedf 


dfdg^ 


dfduj 


dfdi 


dp 



\ 



(6.18) 



taken at the critical points in question (notice that this is not a Hessian matrix, but the 
standard argument for classifying fixed points leads us to handle it as such). After this 
numerical work we find that in all cases where we had AF without gravity we also have 
some solution giving AF with gravity. This is easy to understand by examining the values 
of the fixed points: of all the solutions found for each (A^, ni, n2)-combination, there is at 
least one whose value of ^ tends to be around 1/6 and whose value of / is fairly close to /2 
— the AF fixed point without gravity — which gives AF in the presence of QG. In addition, 
we also find cases where QG makes possible the existence of AF solutions which are banned 
without QG (as was observed in [Q). These new situations correspond to the values of 



16 





rii 


n2 


7 





12 




1 


7 




2 


2 


8 





14 




1 


8 




2 


2 


9 





16 




1 


9 




2 


2 


10 





18 




1 


10 




2 


2 



Table 1: Values of iV, ni and n2. 



N,ni,n2 in the following table and depict a shift in the boundary of the (A^, rii, ^2) -region 
for which AF existed without QG. More precisely, it is a shift by decreasing the allowed 
value of n2 in one unit, as already commented in [Q. Thus, we have shown that the 0{N) 
GUT under discussion, interacting with quantum i?^-gravity, may be considered a completely 
asymptotically free theory (for some given field contents of this model). This study has been 
carried out by introducing a potential in the space of RG couplings, similar to a c-function. 

7 Running gravitational coupling in asymptotically free 
0{N) GUT with quantum i?^ gravity 

Now, having at hand the asymptotically free 0{N) GUT interacting with quantum 
gravity ( |6.1| ) we may discuss the behaviour of the gravitational coupling constant. Unlike for 
GUTs on classical gravitational backgrounds, we cannot analytically solve the RG equations 
for the gravitational coupling constant. Instead, we have a system of RG equations for the 
massive couplings m^, 7 = l/(167rG') and A, which may be analysed only numerically (after 
the corresponding study has been done for the massless coupling constant). This system 
may be explicitly written using the calculations of one-loop counterterms for massive scalars 
interacting with gravity in Ref. and the calculation of the scalar 7-function in Ref. 
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3[. For the Lagrangian (|6.1| ) we obtain 

= /3a = 



dA 



m 

dm? 
~dF 



/3, 



K 



-K 



K 



AA /43 20 1 

\ u H 

2 V 3 3 6cj 

, /lO 1 13 
A7 UJ 

' V 3 6 

rn" (^^— /--(iV-iy 

,2/43 5 ^ 9.3. 

^7 411; 2 



iVm^ ( e - ^ 



(7.1) 



6 



12cj 



1 



33 
2^^ 



7.1 Without gravity 

Ignoring all the QG pieces, and performing the changes (|6.2|) , we are posed with a system 
of differential equations consisting of ( |6l5|) plus 



dA 

dr 

d'j 

dr 

dm? 

~~d^ 



N m . m I 
-— 7 + AiV— e 
4 g 



(7.2) 



It is possible to estimate the type of 7 solution to this system in the large-r regime. 
Considering the most interesting case i.e. A > 0, — and therefore /(r) given by ( |6.14|) — we 
take approximations of the type e^"^ — 1 ~ e^"^ and, using ( |6.11| ), arrive at 



m^ir) 



7(r) 



^(0)e-^^/2, 
7(0) +71 (l-e(^'-m) 



T 



B 



7i 



A^ + 2 



lb^-S(^N-l)]+3iN-l), 



iV(e-l) m^(0) 



_ (7.3) 

Actually, B is positive whenever we are in the settings of AF for both g"^ and / — i.e. 
(A^, rii, 77-2) combinations in the first table — while 6^ — B/2 is negative in these same cases. 
As we see, under the present assumptions 7(r) would tend asymptotically to a constant 
value of 7(0) + 7i as r ^ 00. Undoing now the variable change (|6.2| ), we may put r = 
p- log(l + nb"^ g"^ {0)t) , and writing G(t) in terms of ■y{t) we are left with 



G{t) = G{0) {1 - 167rG(0)7i [(1 + i^b^g\0)ty^l 



r 



(7.4) 
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where, obviously, G{0) = l/(167r7(0)). This expression coincides with ( p.8|) after making 
the notational replacements Ng ^ N, B"^ ^ b"^, 2b + 1 ^ 1 — The above remarked 
asymptotic behaviour is clearly appreciated in Fig. 1. 



7.2 With Gravity 

We can solve numerically the whole system of differential equations and examine the asymp- 
totically free regime of its solutions. By plotting the gravitational running coupling, we 
obtain the curves a and b in Fig. 2. 

8 Gravitational field equations with GUT quantum cor- 
rections in curved spacetime 

Let's now turn to some other application of the running coupling constants in curved space- 
time, namely to the RG-improved effective Lagrangian ( p.3|) . We will work for simplicity 
on the purely gravitational (almost constant) background, supposing that for all quantum 
fields we have zero background. In actual one- loop calculations, it turns out that, working 
in configuration space, the RG parameter t is typically of the form 

1 CiR + 02771^ 

t-^log^-^, (8.1) 

where Ci, C2 are some numerical constants, and is the effective mass of the theory. In 
the model under discussion we have a few different masses, so there is no unique way of 



choosing only one functional form for t [12, 44]. Hence, we will consider the regime of 



strong curvature when curvature is dominant in (|8.1| ). Then, the natural choice of t in the 
RG-improved Lagrangian is (see also |T^, ^) 

*=5log^. (8.2) 

Note that such a regime may lead to curvature- induced asymptotic freedom ||^. In this regime 
we obtain that expression ( |2.3| ) gives the leading-log approach to the whole perturbation 
series. 

= / d^xJ^ I A(t) 

V y I V ; i67rG(t) (8.3) 

+ai{t)R^ + a2(t)C^../3C^'^"^ + a^{t)G] . 

We will consider only asymptotically free GUTs, where A(t), G{t) are given in Sect. 2, 
t = |log^, a2,3(t) = 02,3 + 0.2,3^ (scc (p.7|)). Note, however, that equation (|8.3|) is very 



19 



general. In particular, it has a similar form in the asymptotically free GUTs with i^^.g^-^vity 
of Sects. 6,7 in strong curvature regime (however, the t-dependence of the running couplings 
A(t), . . . , a^it) is of course different). Notice that the technique developed in Ref. |^ may 
be very useful for the explicit solution of quantum corrected gravitational field equations for 
non-constant curvature. 

Working in constant curvature space 



R 



^ R = — 



(3^ ' /52 

we may rewrite (p.3|) expanding the coupling constants up to linear t terms as follows 



U) 



RG 



const X (3^ <^ A + At - 



1 



t 



2^°^ /5V' 



167r 
16 

J' 



G 



+ (ai + ait)— + (as + a^t)^ 



.5) 



where the explicit form of the constants A, G, ai, 03 is evident from (|2.7] ). 

Now one can write the field equation which fixes (3'^ in terms of the theory parameters: 

fr A . ... . , J l \ G n 



8(3^ 



const X 



- + 2(A + At) 



1 



0-^ 



- ( 16ai + -as 



^.6) 



2/52 



0. 



One can see that there exists a real root of (p.6|) for most asymptotically free GUTs (for a 
massless theory, see also discussion in Ref. 0). Hence, we have got a non-singular universe 
with a metric of the form 



ds'^ = o?{r]) drf — 



dx^ + dy"^ + dz^ 
1 + ^/32 



(8.7) 



where k 



1, 0, —1 for a closed, asymptotically open and open universe, respectively, and 



a(r/) 



'2^3/?; 



+ t£ 



\fkri \ 
2V3f3) 



as the solution of the gravitational field equations with quantum GUT corrections. The 
GUTs quantum corrections are here present in the form of (3, which is defined by GUT 
parameters from ( |8.6| ). In terms of physical time T, the scale factor may be put as 

VSf3 r 



a(T) 



+ ke 



-r/(V3/3) 



(8.9) 



20 



Note that such a non-singular inflationary type (De Sitter) solution in gravity theories with 
higher derivative terms (induced by quantum matter) has been discussed in Refs. ||, ^ |41 



in different contexts (for a very recent discussion on non-singular cosmologies in higher- 
derivative theories, see ||4^). 

Let us consider the closed universe {k = 1) in (|8.7|) . In this case, for some choices of the 



theory parameters one can find an imaginary f3 as the solution of Eq. (|8.6|). Similarly to 



what happens in papers such a solution may be interpreted as a Lorentzian wormhole 
which connects two De Sitter universes. However, this wormhole, that for |/?^| ~ Lp^ has a 
'mouth' of the Planck size, results from quantum GUT corrections. From another viewpoint, 
working with (|8.9|) at k = 1 and imposing the initial conditions a(0) = Rq, d(0) = we get 

a{T) = V3(3cosh(^-^y (8.10) 



Some analysis of this solution is given in |^5|. For example, when the effective cosmological 
constant is zero, the solution ( ^.10|) corresponds to a closed universe connected through 
a wormhole to fiat space. In a similar fashion one can construct the quantum corrected 
gravitational equations in other regimes and study their solutions. 



9 Discussion 

In the present paper we have discussed the running of the gravitational coupling constant 
in asymptotically free GUTs in curved spacetime, in the effective theory for the conformal 
factor and in asymptotically free i?^-gravity interacting with an 0{N) GUT. The running 
gravitational coupling constant has been used to calculate the leading quantum corrections 
to the Newtonian potential. These corrections have logarithmic form in asymptotically free 
GUTs and in the effective theory for the conformal factor, and power-like form (but of 
different nature from that in Einsteinian gravity) in finite GUTs. In i?^-gravity with 0{N) 
GUT, the behaviour of the gravitational coupling constant is numerically analysed. Its decay 
rate gets higher as the value of m^(0) is raised. 

The running coupling constants are also necessary in other respects, particularly in the 
RG-improvement procedure. We have found the RG-improved effective gravitational La- 
grangian in the regime of strong constant curvature, and have discussed the non-singular De 
Sitter solution of the corresponding quantum corrected gravitational field equations. The 
present technique is quite general and can be applied in various situations, in particular 
for the construction of RG-improved non-local gravitational Lagrangian, what we plan to 
discuss elsewhere. 
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The other interesting field where the running gravitational constants calculated in this 
paper play an important role is in the quantum corrections to the Hawking-Bekenstein black 
hole entropy. The black hole entropy (Bekenstein-Hawking formula) has the following form 

s-A 

where A is the surface area of the event horizon and G the 

gravitational constant. It has been suggested by Susskind et al. that in the brick 
wall approach (e.g. the 't Hooft cut-off regularization |47]) the quantum corrections to the 



black-hole entropy can be absorbed in the above formula as a simple renormalization of the 
gravitational constant G. In this respect, taking into account the results of our study we 
are led to conjecture that in the case of dimensional regularization (where only logaritmic 
divergences are important) the quantum corrections to the Beckenstein-Hawking formula are 
precisely given by the standard renormalization of G, as discussed above in the present paper 
— plus some less essential contributions from the higher-derivative terms in (2.1). This gives 
further relevance to our calculations in the paper, since then, in order to take into account 
GUT contributions to the black hole entropy one just has to use G{t) instead of G in the 
formula for the entropy. 

Note that the quadratic terms in the lagrangian contribute the black hole entropy already 
at tree level as shown in Refs H3, M. For the charged black hole they give corrections 



hke 

^ — TTT ~i A r ••• 

AG Ah 

From the results of our paper we can futher infer that the one-loop corrections will have the 
form of this equation but with all the coupling constants being now a function of t. We hope 
to return to this question in near future. 
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Figure captions 

Fig. 1. Running gravitational coupling G(t) = l/{16TT'~f{t)) obtained by numerical integra- 
tion of the full system made of Eqs. ( |6.3| ) (written in terms of the original RG parameter t 
of Ref. |3^) and Eqs. (|7. 1|) , for N = 7, rii = 0, n2 = 13. The initial values are g^{0) = 0.1, 
uj{0) = 0.1, ^(0) = 0, /(O) = 0.5, A(0) = 0.1 ,7(0) = 0.1, m^{0) = 0.1. and A(0) = 0, i.e. 
initially there is no quantum gravity. The asymptotic tendency of G{t) towards a constant 
value is clear. 

Fig. 2. G{t) in presence of QG with the same initial conditions as in Fig. 1, except for: 
a) A(0) = 0.5 and m^(0) = 0, b) A(0) = 0.5 and m^(0) = 0.5. The running gravitational 
constant is quickly decreasing. 
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